J/?/) suppression in heavy ion collisions by quark momentum diffusion 
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The momentum diffusion effect of the quark pair due to the multiple scattering in a nuclear medium 
is studied to explain the observed i/ij) yields in SPS experiments. The resulting suppression is found 
to be insufficient to reproduce the J/'i/> yield in Pb-Pb collisions at SPS energy. 
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The suppression of the J/?/' cross section in ultra- 
relativistic nucleus- nucleus collisions has drawn much 
attention since its proposal as a possible signal of the 
quark-gluon plasma (QGP) [1], and the interest is even 
more increased after the anomalous suppression observed 
in Pb-Pb interactions at SPS energy [2]. There have 
been so many works done for explaining this anomalous 
suppression in the conventional manner or in the QGP- 
presuming way [3,4]. 

Qiu, Vary and Zhang (QVZ) recently proposed a 
model [5] for calculating the J/?/; suppression caused by 
the nuclear effect in high-energy proton-nucleus (pA) and 
nucleus-nucleus {AB) collisions. Their model is based on 
the QCD factorization of the inclusive processes. The 
nuclear medium effect is taken into account as a modi- 
fication of the transition probability of the cc pair due 
to the multiple scattering in the nuclear environment. 
The multiple scattering effects have been observed as the 
momentum imbalance of two jets in hadron-nucleus colli- 
sions [6,7]. For the i/ip production in heavy ion collisions, 
the multiple scattering will increase the relative momen- 
tum of the pair and accordingly push the invariant mass 
toward the open charm threshold as it traverses the nu- 
clear medium. QVZ model successfully explains the ob- 
served suppression from pA to AB collisions including 
the Pb-Pb data as nuclear effect [5,8]. 

In QVZ model, the net effect of the multiple scatter- 
ing is expressed as a shift of the relative momentum of 
the quark pair in the transition probability function(see 
Eq. (3)). This shift operator may be obtained by sum- 
ming up all twist contributions at the leading order in 
the strong coupling constant as and the nuclear size, in 
the perturbative QCD calculation [9,10]. This is quite 
different from the result of the classical multiple scat- 
tering process; after independent multiple collisions with 
random momentum transfers, the distribution of the rel- 
ative momentum q of the pair diffuses while the mean 
increases. The only aim of this Note is to examine 
to what extent the J/V' yield in heavy ion collisions is 
explained as the influence of the momentum diffusion, 
treating the multiple scattering of the quark pair as ran- 
dom walk. We show that the resulting suppression is 
insufficient to explain the Pb-Pb data and looks rather 
similar to the Glauber model result. 

Let us first review QVZ model [5] . The inclusive cross 
section of the J/'0 production in the collision of the 



hadrons A and B is written in a factorized form: 
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i^cc->J/^(?'), (1) 



where t runs over all parton fiavors, = q^ + 4m^, 
(j)a/Ai^a) is the distribution function of parton a in 
hadron A, and xp — Xa ^ Xb and XaXb = Q^/s. The 
parton cross section a is given in [10]. This is the lead- 
ing order formula in and the phenomenological con- 
stant Kj/^ corrects the higher order effects. F^^^j/^{q^) 
describes the transition probability for the cc state of 
the relative momentum q^ to evolve into a physical J/i/i 
meson. For this transition probability they propose a 
parametrization of 



^^i2j/j9') = N„^e{q^)6{Am'^ - Ami ^ 



4m' 



4to? 



(2) 



which includes the effect of the open charm threshold 
at 4m'^ and simulates the gluon radiation effect with 
the parameter > with putting the larger weight 
to the smaller . The constant probability [ap = 0) 
corresponds to the color-evaporation (CE) model. 

For the J/^ production in the y>A and AB collisions 
QVZ model additionally assumes the separation of the 
multiple scattering which affects the pair state and the 
formation of the J/V' resonance in the high-energy in- 
teractions. The multiple scattering of the pair in the 
nuclear medium would increase the relative momentum 
q^ of the pair. The effect of coherent multiple scattering 
in the perturbative QCD calculation may be represented 
by shifting of the relative momentum in the transition 
probability [5] as 



(3) 



where L is the effective length of the nuclear medium in 
the AB collisions. We note here that for a large enough 
L such that (f > 4m'^ — 4to^ the transition probabil- 
ity essentially vanishes due to the existence of the open 
charm threshold (2). This apparently gives rise to a much 



1 



stronger suppression than the exponential one following 
from the Glauber model. 

We would like to describe now the momentum diffusion 
effect by modifying the transition probability i^cc^J/i/- 
of QVZ model. The cc pair with relative momentum, 
q, produced in a hard parton collision, will change its 
momentum to q' after the random multiple scattering, 
and then transforms into the J/V' with the probability 
Fcc^}/^{q''^). Here we treat the exchanged soft momenta 
as independent and random in the multiple scatterings. 
After many scatterings, this classical, elementary diffu- 
sion process of the momentum results in the Gaussian 
distribution around the initial value q with the variance 
e^L. The three-dimensional random walk is assumed here 
only for simplicity. The transition probability should be 
in effect replaced by the one smeared with this Gaussian 
weight as 

(4) 

Note that the transition probability i^cc^j/j/'('?^) never 
vanishes for any q although the average momentum of 
the pair increases as {q'^) = q^ + 3e^L after the multi- 
ple scattering. We immediately find that the transition 
probability behaves for the asymptotically large L as 

Fc,^j/Al') - -(2^J^3j2 I d\'F,,_^„^{q'\ (5) 

This power suppression in L which stems from the deple- 
tion of the normalization factor is more moderate than 
the exponential one. 

In Fig. 1 we show our result on the J/'^ suppression 
calculated using the formula (1) with the smeared prob- 
ability (4) . The parton distribution function of CTEQ5L 
[11] is used without any nuclear modification, and the 
parameters are fixed to the same as [5]: Q!_f = 1 and 
h/4' = K:]i^,N:]i^, = 0.485 for the transition probability 
(2). The multiple scattering is included as smearing with 

= 0.185 GeV^/fm. Our model reasonably fits the data 
in the p^ and AB collisions taken from [2] except the Pb- 
Pb point. The curve bends upward in the semi- log plot 
as is expected from Eq. (5), but is almost consistent with 
the straight line within this interval of L. The original 
QVZ model (3) with = o.25 GeVVfm (dashed line) 
can explain all the data points in Fig. 1. The downward 
bending of QVZ model is the result of the existence of 
the open charm threshold in the transition probability (2) 
and the uniform shift of the momentum (3). 

We performed the calculations using other forms for 
the transition probability, the Gaussian form and the CE 
form, besides Eq. (2) [5], and confirmed that the quali- 
tative behavior of the suppression is unaltered; the Pb- 
Pb data alone lies far below the curve calculated by our 
model. The concavity of the suppression curve in the 
semi-log plot is a robust consequence of the momentum 
diffusion Eq. (4) by the multiple scattering, irrespective 




2 4 6 8 10 
L(A,B) [fm] 

FIG. 1. 3/ip cross section to n'^ fJ,^ pair [2] compared with 
the models using the smeared transition probability Eq. (4) 
(solid line) and the transition probability with the shifted 
momentum Eq. (3) (dashed line). 

of the detailed form of the transition probability. There- 
fore the Pb-Pb data cannot be explained as the momen- 
tum diffusion effect of the quark pair due to the multiple 
scattering before forming a physical resonance. 

We like to comment that in the high-energy collisions 
the transverse momentum diffusion in two dimensions 
might be more appropriate than the three-dimensional 
one, as in the Glauber approach. Then the L dependence 
of the suppression would change from to L~^. 

The absorption of a hadronic state by a power in L 
is similar to the result of the color transparency at high 
energies predicted in many theoretical studies [12,13]. In 
the model studied here, however, the color degrees of 
freedom of the pair states are not explicitly treated in the 
transition probability nor any interference effect. Hence 
the relation to the color transparency effect is unclear 
within this model. At the collider energies like RHIC 
and LHC the coherence effect becomes more important 
and the quantum description of the resonance production 
at parton level should be more elaborated. 

The uniform shift of the relative momentum is the re- 
sult of the selective sum of higher twists at the leading 
order of as and the nuclear size [9] . From the view of our 
model, it may be worthwhile to investigate the smear- 
ing effect due to the correction terms in the perturbative 
QCD calculation. 

In conclusion, the suppression of the J/tp cross section 
in AB collisions except the Pb-Pb at SPS energy can 
be described by the momentum diffusion by the multiple 
scattering of the cc pair before the resonance formation. 
The suppression of the model is insufficient to explain 
the J/i/j yield observed in Pb-Pb collisions at SPS energy, 
which is therefore left anomalous in this treatment. 
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